Manipulating unpaired Majorana fermions in a quantum spin chain 



O 



Abhinav Saket, S. R. Hassan and R. Shankar 
The Institute of Mathematical Sciences, C.I.T. Campus, Chennai 600 113, India 

(Dated: July 26, 2010) 

We analyse an exactly solvable spin-1/2 chain which is a generalised version of Kitaev's honeycomb 
model. We show that every state of the system has a 2^' * fold degeneracy, where A'^ is the number of 
sites. We present analytic solutions for the zero energy modes of the Majorana fermions. Localised, 
unpaired Majorana modes occur even in the bulk of the chain and they are bound to kink (anti-kink) 
Z2 flux configurations. The unpaired Majorana modes can therefore be created and manipiulated 
if the Z2 flux configurations can be controlled. We delineate the regions in parameter space for 
homogenous chains where the zero modes occur. We further show that there is a large parameter 
space for inhomogenous chains where the unpaired modes occur and that their wavefunctions can 
be tuned if the couplings of the model can be tuned. 
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I. INTRODUCTION 

The idea of topological quantum computation as a way 
to incorporate fault tolerence at the hardware level has 
been getting a lot of attention recentlj^. In this scheme, 
qubits are non-abelian anyons and the braiding opera- 
tions on them implement quantum gates. One of the 
simplest class of non-abelian anyons are realised in sys- 
tems with unpaired Majorana fermions (UMF)^. In a 
fermionic system with N zero energy modes, there are 
2N Majorana modes. If these 2N modes can be inde- 
pendently moved around each other, then the geometric 
phases picked up correspond to a non-abelian represen- 
tation of the braid group. Thus it is of interest to study 
physically realisable systems where UMF exist and can 
be manipulated. 

Non-abelian anyons are theoretically predicted to oc- 
cur in certain fractional quantum Hall states like v — 
5/32.. There is also theoretical work showing how they 
could be realised in quantum circuits^. In this con- 
text, Kitaev presented a remarkable solvable spin-1/2 
model on a honeycomb lattice^ which realises non-abelian 
anyons made up of UMF. The model can be written 
in terms of Majorana fermions in the background of Z2 
gauge fields. This fermionisation procedure is very sim- 
ilar to that proposed in the context of the resonating 
valence bond (RVB) theory-. The crucial difference is 
that in Kitaev's honeycomb model, the gauge fields are 
constants of motion, hence the RVB type mean field the- 
ory is exact and the problem reduces to solving a theory 
of non-interacting Majorana fermions in the background 
of static Z2 gauge field configurations. Kitaev showed 
that the ground state is the flux free configuration. The 
model has a phase which is characterised by a topologi- 
cal invariant, the Chern number, being equal to ±1. In 
this phase there are UMF trapped to each vortex (a Z2 
flux). Thus if vortices can be created and manipulated, 
it is possible to braid the UMF. However, the Z2 flux 
operator is a 6-spin operator and thus may not be easy 
to realise in practice. 

Kitaev's honeycomb model can be generalised to a va- 
riety of other lattices. It can be constructed on any lattice 



with coordination number three, if all the bonds can be 
coloured using three colours. It has been shown that all 
such models can be realised in cold atom systems^ and 
using quantum circuitsS. 

In paper we present a one dimensional generalisation 
of Kitaev's honeycomb model which we call the tetra- 
hedral chain (TC) and analyse its zero energy modes in 
detail. One dimensional models with UMF at the edges 
have been studied earlieriSiiS. The new feature of the 
TC is that the wavefunctions of the UMF not necessarily 
peaked at the edges of the chain but can be peaked any- 
where in the bulk. As we will show, they are trapped to 
kink and anti-kink flux configurations and can be moved 
by tuning the fiux configuration. Further, by tuning the 
coupling constants, their wavefunctions, which we obtain 
analytically, can also be tuned. 

The rest of this paper is organised as follows. In sec- 
tion|lT]we present the model and its conserved quantities. 
Section IIIII describes the Jordan- Wigner transformation 
which enables us to rewrite the theory in terms of Ma- 
jorana fermions hopping in the background of static Z2 
gauge field configurations. The diagonalisation of the 
Majorana fermion problem is described in section ITVl and 
the numerical calculations to determine the ground state 
fiux configurations are presented in section [V] Section 
I VII gives an analytic proof of the degeneracy of the eigen- 
states that we find numerically. The dependence of the 
gap as a function of the couplings is computed in section 
IVIII and a gapless line in the parameter space is iden- 
tified. Section IVIIII contains a detailed analysis of the 
zero-modes of the Majorana fermions. We summarise 
our results and discuss them in the conluding section [ 



II. THE TETRAHEDRAL MODEL 



The chain we define our model on is shown in Fig. [T] 
There are four sites per unit cell which are labelled as 
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FIG. 1; The tetrahedral Kitaev chain. There are four sites 
per unit cell. The x, y and z bonds are as indicated. 



shown in the figure. The hamiltonian is, 

^ = I]('^4<-l,4<l+<2<3) 



'^\,-i'^X^) 
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There are two triangular plaquettes in each unit cell 
and a conserved Z2 flux associated with each of them. 
The flux operators are, 



\mL ^^ ^z ^y T-r/^ ^^ ^z ^y 



(2) 



As in Kitaev's honeycomb model, these quantities are 
conserved as a consequence of a local spin rotation sym- 
metry of the model. Namely, a it rotation on each site of 
a plaquette about the direction of the outgoing bond. 

Apart from these local conserved quantities there are 
also three global qauntities which are conserved as a con- 
sequence the fact that a global tt rotation about each of 
the three axes is a symmetry of the model. We denote 
these by, 



Y." EEe^^^'-<' 



(3) 



It can be verified that E^ is the product of the fluxes of 
aU the plaquettes, 

S- = n {WtW^) (4) 



III. FERMIONISATION 

We express the hamiltonian in terms of Majorana 
fermions using the Jordan- Wigner transformation. We 
choose the Jordan- Wigner path to go along the x and 
the y bonds from left to right. At every site we have two 
bonds that are tangential to the path. We denote the in- 
coming bond by ti . This is the a;-bond on the sublatttices 
1, 3 and y-bond on sublattices 2,4. The outgoing bond, 
corresponding to y on sublattices 1,3 and x on sublattices 
2,4, is denoted by ^2- The third bond on each site which 
in normal to the path is denoted by n with the sign de- 
fined hy fi — tixt2- With our choice of path, the normal 
bond is z for sublattices 1,3 and —z for sublattices 2,4. 



With these definitions the two Majorana fermions at 
each site are defined as. 
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The hamiltonian can be expressed in terms of the Majo- 
rana fermions defined above, 

H ^ ^{Jx {iii-iA^i,! + «6,26,3) 
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where the link fields, u 
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are defined as, 
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(9) 



It is easy to see that the link fields are conserved quan- 
tities. Thus, as expected for a generalised Kitaev model, 
the theory gets written in terms of Majorana fermions 
with nearest neighbour hopping in the background of 
conserved Z2 gauge fields with the gauge fixing condi- 
tion that the gauge fields on the x and y bonds are equal 
to -f 1. It is easy to check that in this gauge the two flux 
operators are proportional to the link fields. 



"yi,3??j,i = W^ 
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(10) 



It is instructive to express the three global conserved 
quantities, S" in terms of the fermionic variables. We 
have. 
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(11) 

(12) 
(13) 



We will refer to E^ as the flux number and S^ as the 
Majorana number. 

IV. DIAGONALISATION 

The hamiltonian can be diagonalised in the standard 
way. We write the eigenstates as direct products of states 
in the 77 fermion sector, \Q)^ which we will refer to as the 
gauge sector and states in the ^ fermion sector, \M), 
which we call the matter sector. We choose the states 
in the gauge sector to be the simultaneous eigenstates of 
the Z2 flux operators, i.e \Q) = |{wf',wf-}), where 

«f ^""^K^f ,^f }) = ^f^^'^K^f ,^f }) (14) 



We then have 

H [u^.uf] \M)\{wtw^}) = H [wtwf] \M)\{wtwf}) 

(15) 
The problem reduces to finding the eigenstates of the 
quadratic hamiltonan of the ^ fermions in the background 
of the gauge field configuration {wf ,wf }. The normal 
modes are given by the solution of the eigenvalue equa- 
tion, 



3P 



(16) 



where A = T + y is a purely imaginary anti-symmetric 
matrix, 
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(17) 



(18) 



The eigenvalues come in pairs, ±e„. So we can always 
choose e„ to be positive. The eigenvectors correspond- 
ing to the positive and negative eigenvalues are complex 
conjugates. Denoting the real and imaginary parts of the 



eigenvectors as 



inBil) 



we expand the ^ fermions as. 



6:.a = Y. """^^^ + ^ 



'nVi,a 



(19) 



a„ and &„ are also Majorana fermions. The diagonal form 
of the hamiltonian is then. 



H 



E 



(20) 



The ground state energy is the one where all the oper- 
ators ibnttn are diagonal and equal to -1. The ground 
state energy is. 



Eo 



E 



(21) 



A. Boundary Conditions 

We will be analysing the system with periodic (PBC) 
and open boundary conditions (OBC). The latter case is 
straightforward. The fermionic hamiltonian is exactly of 
the form given in equation ([7]) with i = 1, .., iV, where N 
is the number of unit cells. The fermionic modes have to 
be solved with the boundary conditions 



i>0,a = 4>N+l,a — 



(22) 



PBC is a little more subtle. As is standard in Jordan- 
Wigner transformations, the term in the hamiltonian for 
the link i = A^ to i = 1 is. 



^JV,0 — Jx^^i^l,l^N,4 



(23) 



S^ is a conserved quantity and can hence be chosen to be 
diagonal. The fermionic modes have to be solved with pe- 
riodic boundary conditions for states with S^ = -1-1 and 
with anti-periodic boundary conditions for states with 
S^ = -1. Namely, 



t>N+l,a 



(24) 



where p — ±1 is the eigenvalue of S^ 



V. THE NUMERICAL SOLUTION 

We have solved eigenvalue problem numerically for a 
chain with N unit cells and open boundary conditions for 
N < 5. Each unit cell has 4 possible flux configurations 
making a total of 4^. We computed the single particle 
spectrum for each flux configuration and calculated the 
lowest total energy state for the range of values Jx = 
1.0, Jy = 0.1 - 10, Jz = 0.1 - 10 (in the steps of 0.2). 

We find that the ground state is 2^ fold degenerate 
for all values of the parameters. This denegeracy comes 
from the fact that the single particle spectrum depends 
only on the values of the product, Wi = wf"wf, and not 
on their individual values, i.e the spectrum only sees 
the total fiux passing through each unit cell which is the 
product of the values of the fiuxes passing through the 
two triangular plaquettes that make up the unit cell. 

The ground states correspond to the translationally 
invariant, Wi = 1, flux configurations. 



VI. THE DEGENERACY OF THE STATES 

As mentioned above, we find that the spectrum in each 
sector is 2^ fold degenerate. This degeneracy is due to 
the fact that the energy eigenvalues depend only on the 
total flux in the unit cell, namely Wi = wf'wf' and not 
on the individual fluxes, wf" and wf. In the remaining 
part of this section, we will give an analytic proof of this 
2^ fold degeneracy of all the states. 

This degeneracy is related to but not the same as the 
4^ degeneracy that occurs in the simple Kitaev chain, the 
Jz = limit of our model. In this case the degeneracy 
is easy to understand. The gauge fields do not occur at 
all in the hamiltonian. Thus each state is 4^ degenerate 
corresponding to all the states in the gauge sector. The 
extra z-bond terms in our model lift this degeneracy only 
partially. 

At Jy — Jz, the denegeracy can be understood in terms 
of a local symmetry. It consists of interchanging the spins 
at sublattice 2 and 3 in any unit cell and then performing 
a 7r/2 rotation about the x-axis on all the spins in that 



unit cell. The operator that implements this transforma- 
tion is, 



R^(^^i:^^i±l]e^f5:t^i< 



(25) 



It is easy to verify that the Pi 's commute with the hamil- 
tonian, However Pi flips the sign of the two flux operators 
in the unit cell, 



p,W.^(«)i^, 



-W,^(^\ Pf^Wtwl" (26) 



Thus it changes the flux conflguration while conserving 
the total flux through the unit cell . Since it also does 
not change the energy eigenvalue, it follows that every 
eigenstate of the hamiltonian is 2^ fold degenerate. 

However, we numerically observe that the degeneracy 
persists even when Jy ^ Jz- We will now give a proof 
for the degeneracy which is valid at all couplings. We 
note that in equation P^ . (j>i^2 and (j)i^3 couple only to 
sites within the unit cell. We express them in terms of 
4>i^i and (f)i^4 and obtain an eigenvalue equation for these 
quantities. We are then able to show that the eigenvalues 
depend only on Wi. 

We deflne the two component column vectors. 
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and the matrices. 
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The eigenvalue equations (J16p can be written as. 
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(27) 

(28) 
(29) 

(30) 



where r", a = 1,2,3 are the Pauli matrices. ip can be 
eliminated from the equations to get. 



T + U- 



1 



rC/t 



X = eX 



(31) 



e - JxT^ 
The equations can be explicitly written as, 

- iJ2,(/)i_i,4 + Cie''"'(/>i^4 = X(f>i^i (32) 

iJx4'i+i,i+ Cie~''°''(l)i^i = \4>iA (33) 

where. 



When these equations are solved for A, they will yield an 
equation for e. We now make a transformation. 
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This gets rid of the phases in equations ([5^ and ([55)) 
which become, 

- iJx4>i~iA ~ 'ici(f>iA = A0i,i (40) 

iJx(t>i+i,i + ici(f>is = A0i^4 (41) 

Since Ci depends only on Wi and not on w^ and wf' indi- 
vidually, A and hence e depends only on the Wi . Thus the 
energy eigenvalues depend only on the total flux passing 
through the unit cell. 

This result is true for all values of Jx , Jy and Jz . Note 
that when Jz = 0, equation (j34l) implies that Ci is inde- 
pendent of Wi also. Thus the 4^ fold degeneracy of the 
simple Kitaev chain is recovered in this limit. 



VII. THE GROUND STATES AND GAPS 

As we have mentioned earlier our numerical results 
show that the translationally invariant fluxes through the 
unit cells, namely Wi = 1 is one of the ground states of 
the model. The fermionic problem is easy to solve an- 
alytically for these flux configurations and therefore the 
energy of the w = 1 configuration can be compared to 
that of the w = — 1 configuration in the thermodynamic 
limit. 

If u;i = w, then a in equations (j40l[4T|) are independent 
of i, Ci = c. The equations can be solved by Fourier 
transforms. 






e-*>..a 



(42) 
(43) 



A is then given by the eigenvalues of the following matrix, 

(44) 
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namely, 
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(35) 
(36) 



A = ±y/c'^ + Jl +2cJxCosk 



(45) 



Along with equation p6p this yields an equation for the 
four energy bands. We have solved for this energy bands 
and have computed the gap for < Jy,Jz < 10 and 

Jx^l- 

The calculations confirm that the Wi = +1 sector has a 
lower ground state energy than the Wi = —1 sector. The 
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FIG. 2: The gap plotted against Jy for Jo; = 1 and values of 
Jz varying in the range < Jz < 10. It can be seen that the 
Jy < 1 region is gapless for some value of Jz. 



fermionic gap is twice the value of the lowest single par- 
ticle energy eigenvalue. This is shown in Figure [2] where 
the gap is plotted against Jy for different values of Jz- 
The Jy < 1 region is clearly gapless. It is interesting that 



when we plot the gap as a function of J = , / J^ + J|, all 

the points fall on the straight line as shown in Fig. |31 
This can be shown analytically since the expression in 
equation (|34|) simplifies considerably at Wi = 1. It is 
then easy to show that the energy eigenvalues are given 
by, 



J^ + l±2Jcos( - 

The gap is thus given by, 

A = 2JJJ-1| 



(46) 



(47) 



The circle in the parameter space, J^ ^ Jy + J^ is there- 
fore gapless. 



VIII. ZERO MODES AND KINKS 

We will now analyse the zero energy modes of the Ma- 
jor ana fermions. When e = 0, we have A = and hence 
the equations ([iU)) and dH]) decouple and become very 
simple, 



Jx4>i+l,l + 
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FIG. 3: The gap plotted against J = i/J| + J|. The points 
fall on straight lines with slopes ±2 



These recursion relations can be formally solved. 
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Where (/ii(4) are arbitrary constants. Thus, there are two 
formal solutions for every set of values of the parameters 
and every flux configuration. One with (\)i^\ ^ and 
(/)i,4 = which we denote by (j)^ and the other with (pi^i — 
and (jji^A ^ which we denote by (j)~ . 

However, the boundary conditions that the modes have 
to satisfy will pick out certain flux configurations for each 
point in the parameter space. We will analyse the situa- 
tion for the cases of periodic boundary conditions (PBC) 
and open boundary conditions (OBC). 



A. Periodic Boundary Conditions 

We consider a chain with N unit cells. As discussed in 
equation ([24]), PBC will imply that. 



4'N+l,a =P4>l,a 

Equations 0, dSB) and dH]) imply 
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Consider the general case where AI < N of the w^'s are 
equal to —1 and iV — M of them are -1-1. We will refer 
to such configurations as M-defect configurations. Equa- 
tion (|53p gets written as, 
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(54) 



When Af = 0, we have only have solutions on the circle 
of radius J^ in the Jy — J^ plane. Note that M = is the 
translationally invariant ground state flux configuration. 
Thus this result implies that the model is gapless only on 
the circle, consistent with equation (|47l) . 

When A/ > 0, since the LHS of equation ((54|) is less 
than 1, no zero modes exist within the circle of radius 
Jx- Outside this circle, for every even M there are two 
directions where equation ([54]) is satisfied with p = +1. 



One direction where Jy > Jz and the other where Jz > 
Jy. For odd M, the Jy > Jz solution exists for p = +1 
and the Jy < Jz solution ior p — —1. 

Thus for every N there are a discrete set of points 
outside the circle which support zero energy modes. In 
the thermodynamic limit of iV -^ oo, N/M can take all 
values from 1 to cxd. In this limit, all the points outside 
the circle in the range, J^ < \Jy — Jz\ < support zero 
energy modes. This region is shown for Jj; = 1 in Fig. 2] 



6j 2 a-nd (f>^ 3 are given in terms of (j)^ -^ and 0^ ^, 
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It can be seen that 0+ is peaked ai i — M and is 
minimum at i = 1 whereas 0~ is peaked at i = 1 and 
has a minmum at i = M. Thus we have one Majorana 
mode localised at the location of the kink and another at 
the location of the antikink. When M is large these are 
well separated. Thus the situation is similar to Kitaev's 
honeycomb model with kinks and antikinks playing the 
role of the vortices. If the flux configuration can be ma- 
nipulated, then so can the Majorana modes trapped to 
them. 




FIG. 4: The region in the Jy — Jz plane at J^ = 1 that 
supports zero energy modes 

Consider the cases when the M defects are in adjoining 
unit cells, say from i = 1 to z = M. We will call this a 
kink-antikink configuration. We define polar coordinates 
in the Jy — Jz plane. 



■^ = \ Jy + -^z, 7 = tan" 
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The unnormalised wave functions of the two zero modes 
are, 



't. 
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= j2(*-i) i<M + l 



(56) 



J2(i - 1) (cos27)*+'~*'^ i>M+l (57) 
(58) 

j-2{t-i) i<M + 1 (59) 

J-2('-i) (cos 27)-('+i-*') t>M + l (60) 
(61) 



B. Open Boundary Conditions 
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FIG. 5: The wavefunctions of the two Majorana zero modes 
for A*' — 25, M = 7 with open boundary conditions. 



We now consider open chains with TV unit cells. We 
then need to solve equations (|48| and (|49l) with the 
boundary conditions in equation (j22[) . From the solu- 
tions in equation ((50|) and (|5ip . we see that the above 
boundary conditions have non trivial solutions if and 
only if at least one of the factors in the products on the 
RHS of the equations is zero. This is only possible when 
\Jy\ ~ \Jz\- Thus zero modes exist in this case only at 
7 = 7r/4. At these points, the zero modes are similar 
the ones in the PBC case except that the wavefunctions 
strictly vanish in the region between the kink and the an- 
tikink. Namely, if we consider a flux conflguration with 
Wi — —1, ii < i < ii -\- Af , then the wavefunctions are 



given by, 

<i = J'(^"') ^<h (64) 

= i>ii (65) 

<4 = (66) 

(/)-4 = J-2(^-i) z>zi+M (67) 

= i> M+1 (68) 

'/'M = (69) 

(f>^ 2 ^-iid '/'j 3 ^^^ given in terms of (j)^ -^ and 0^ ^ as before. 
These wavefunctions are shown in Fig. [5] for A^ = 25 and 

M = 7. 

C. Inhomogenous chains 



The solutions for the zero mode equations (j48)) and 
(|^5)l hold even for the case of inhomogenous chains where 
the coupling constants Jy and J^ depend on z. The equa- 
tions then read 



^z ' ^i zi 

Jx(j)i-IA ^ j (j)iA = 

^x 



(70) 
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These solution to these recursion relations is exactly the 
same as the case of the homogenous chain with i depen- 



dent J, 
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(72) 
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It is clear that by tuning the values of the site dependent 
couplings, a large variety of zero mode wavefunctions can 
be engineered. 



D. Tuning the flux configuration 

We can add the following "chemical potential" term 
for the conserved flux operators, 



^. = E (Mf WA^ + Mf W,«) 



(74) 



then any particular flux configurations can be made the 
ground state. 

However it is still not known how to engineer these 3- 
spin operators in the physical realisations of the model 
in cold atom systems or in quantum circuits. 
IX. CONCLUSIONS 

To summarize, we have analysed an exactly solvable 
spin-1/2 chain which we call a tetrahedral chain. This 
model is a generalization of Kitaev's honeycomb lattice 
model. Like the honeycomb model, this model too has 
three coupling constants, J^, Jy and J^. One of them 
can be scaled away, so without loss of generality, we can 
set Jj. = 1. 

The model has conserved Z2 fluxes on every triangu- 
lar plaquette, namely two fluxes per unit cell. We de- 
note the values of these fluxes hy w^ . The fermionic 
spectrum in the background of the flux configuration de- 
pends only on the total flux in the unit cell, wi = w^wf" . 
Consequently, every eigenstate of the model is 2^ fold 
degenerate, where N is the number of unit cells. 

The ground state flux configurations of the model are 
the ones with wi = 1. The model is gapless on the unit 
circle in the Jy — Jx plane and is gapped everywhere else. 

We have found simple analytic solutions for zero en- 
ergy wavefunctions for the Majorana fermions. We have 
delineated the region in parameter space where these so- 
lutions exist. In the uniform case, the two zero modes 
are peaked at the kink and anti-kink positions and are 
well separated if the kink and antikink are well separated. 
The analytic solutions also show how the zero mode wave- 
functions can be engineered by tuning the couplings Jy 
and Jz in a site dependent way. 

In conclusion our results show how zero energy Ma- 
jorana modes can be created and manipulated in the 
tetrahedral spin chain if it is possible to tune the flux 
configurations and the local couplings. In this work, we 
have not addressed the interesting and important ques- 
tion as to how this is useful to braid unpaired Majorana 
modes. While this is of course not a meaningful opera- 
tion in a chain, it may be possible to implement braiding 
in a network of coupled chains using the results obtained 
in this paper. We are pursuing this line of work and hope 
to report on it in future. 
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